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Abstract 

In order to analyze the effect of chaos or order on the rate of decoherence in a subsystem, 
we aim to distinguish effects of the two types of dynamics by choosing initial states as random 
product states from two factor spaces representing two subsystems. We introduce a random 
matrix model that permits to vary the coupling strength between the subsystems. The case 
of strong coupling is analyzed in detail, and we find no significant differences except for very 
low-dimensional spaces. 

PACS: 05.45.Mt, 03.65.Yz 

1 Introduction 

The discussion of decoherence phenomena has been centered around problems of quantum òptics 
and atòmic physics, which usually implied the use of smooth, often Gaussian, wave packets. Yet 
recent ideas on how to process quantum information changed that situation to some extent. Indeed 
in this context it is desirable to store a maximal amount of information in product states of qubits, 
such that it can be processed by unitary time evolution. During the process, the relevant states 
are changed in a complicated way Hence, they will most certainly resemble random states much 
more than Gaussian wave packets. As the conservation of coherence is essential for any quantum 
information application, we have a clear interest in understanding decoherence of random states in 
such situations. The Standard tool in this context is the study of "fidelity" which involves a change 
in the Hamiltonian, and a number of interesting results are available Q, § |. 

The fundamental question, how the integrability or chaoticity of the corresponding classical 
system i.e. "quantum chaos" afïects the process of decoherence (see e.g. ||, [|, ||), appears in a 
diffcrcnt light, because the behaviour of random states becomes relevant. We shall focus on this 
aspect, concentrating our interest on the effects of the chosen dynamics. At first glance, the semi- 
classical context implicit in this view may seem odd after the reference to quantum computing, but 
note that in many-body systems some effects of the chaotic dynamics tend to survive even for very 
low lying states where semi-classical arguments certainly do not apply. 

In a recent letter Q we proposed to discuss decoherence in a random matrix model, in order to 
separate effects of the Hamiltonian from those of the choice of the initial pure state. Such a model 
is ideally suited to discuss the behaviour of random states. In this context we do not embrace 
the scheme of fidelity analysis Q , but simply follow decoherence as a function of time for a fixed 
Hamiltonian. We introduce decoherence by performing partial traces over a subsystem instead of 
considering non-unitary time evolution. What we show, could therefore more strictly be considered 
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as a study of the progress of entanglement of a product wave function under an integrable or chaotic 
time evolution. 

In section ^| we present our random matrix model, which allows for a variety of situations both 
with respect to the nature of the subsystems and the coupling strength between them. Note that 
models based on random matrix theory (RMT) have been used to discuss dissipation (see || and 
references therein), but these aspects will not be discussed here. In section || we shall analyze a 
special situation of strong coupling between the two subsystems. Though not particularly relevant 
for quantum computation, it has the advantage of a transparent mathematical structure which 
allows to obtain non-perturbative analytic results. These are compared to numerical calculations 
in section |J. As it turns out, the results are typical for other situations as well. In particular 
for weak coupling, analyzed numerically in the second part of section [|, we are lead to essentially 
the same conclusions as in the strong coupling case. Our approach to the problem with (largely) 
random wave functions actually separates rather nicely the effects of the dynamics from those of 
the functions considered, and in section ^ we shall conclude with some considerations on this point, 
as well as with a comparison of our findings with those of Refs. (|, |). 



2 The random matrix model 

Properties of chaos and integrability of a classical system manifest thcmsclves both in the spectrum 
and in the wave functions of the corresponding quantum system. While the former is invariant 
the latter are basis dependent. This does not mean that the latter are irrelevant in a semi-classical 
context; they certainly reflect the special features of the dynamics such as KAM tori or short periòdic 
orbits. However, while localized wave packets may feel a strong influence on their dynamics, random 
wave functions would be rather insensitive to such localized features. 

To construct our RMT model we start from the Standard quantum chaos conjecture, that the 
classical ensembles |ï(], [nj (e.g. the Gaussian orthogonal ensemble (GOE) for time reversal invari- 
ant systems) describe the universal features of a quantum system, whose classical counterpart is 
chaotic Jl2| . For classically integrable systems we expect a random quantum spectrum if we exclude 
harmònic oscillators |Ï3| . The random spectrum can be combined with the orthogonal invariance 
to give the so called Poisson orthogonal ensemble (POE) JÏJ], extending in this way the classical 
ensembles to the integrable case. This concept provides the ideal tool to describe the evolution of 
random wave functions, which is one of the main purposes of the present work. 

Consider a Hamiltonian consisting of three terms H = ií (0) + À^ (1 ' 2) with ií (0) = + /i (2) , 
where the two terms of act on different degrees of freedom of the system; we may refer to them 
as the central system and the environment. Note that the foW may act on one or several degrees 
of freedom each, and the interaction í/( 1,2 ) may or may not induce chaos. Indeed the total system 
may be integrable and separable in a different set of coordinates. 

Taking this Hamiltonian as a quantum operator we shall denote by 7íi and TI2 the Hilbert 
spaces of the central system and the environment with dimensions n and m , on which and 
respectively act. The total Hamiltonian H acts on the product space TL = TL\ x TÍ2 with dimension 
N = nm. We write the basis states of Hi as kets with Latin letters such as \i) and those of TL2 as 
kets with Greek letters such as \fx). The states \i, ji) = \i)\fj) with indices conveniently written as 
pairs, form an eigenbasis of H^°\ H is diagonal in a different basis, which we enumerate by a single 
index a. Thus Hí^y l·l > = J2 a ^v*i<* ^ a Oi'v'<a wnere E a denotes elements of the diagonal energy 
matrix E and O the orthogonal transformation between the two bases. 

We distinguish between strong and weak interaction. The interaction strength is usually dis- 
cussed in terms of the spreading width, which indicates the width of the distribution of the expansion 
coefficients of the eigenstates of H in terms of those of H^ 3 ' . This spreading width is a semi-classical 
quantity in the sense that it can be calculated from a phase space integral |ïl| . It is closely related 
to the width of the LDOS (the density of the eigenstates of H expressed in terms of those of H) 
If the spreading width is small, i.e. if the eigenstates of Hq contain a dominant amplitude in the 
eigenbasis of H we clearly have a weak coupling case, where perturbative calculations should yield 
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the correct answer. A detailed study will be published elsewhere, but we shall give some numerical 
results in section [|. 

We now determine appropriate matrix ensembles for the three terms in the Hamiltonian. As 
mentioned before, M 1 ' could describe the central system and the environment, e.g. the heat 
bath. In any case both hP^ and could pertain to either of the above ensembles, the GOE or 
the POE, and y' 1 ' 2 ) could typically be symmetric with independent Gaussian distributed matrix 
elements. 

Finally we must specify the initial states. As mentioned before, we are interested in states of 
complicated structure, but require that they were initially pure. Hence we choose two random states 
in TLi and TI2 and define the initial state as the product state of both. The probability measure for 
the random states is the orthogonally invariant one. 

The situation first studied by Zurek and coworkers |Q] could be simulated as a weak coupling case 
using a GOE with high level density for and a GOE or a POE respectively for However, 
we would have to use initially smooth wave packets. In this respect the states, whose evolution is 
discussed here, are entirely different. 



3 Purity decay in the strong coupling regime 

In the case of strong coupling hP^ and h^ 2 ' determine the factor spaces Tíi and TÍ2 only. Thcir 
spectral properties are irrelevant except for their relative spectral density in the energy region 
where the wave packet lives. The total Hamiltonian which is essentially equal to the interaction 
yí 1 ' 2 ) will be given by the GOE for chaotic systems and by the POE for integrable ones. Both 
ensembles are given by matrices of the form OEO T , where E is a diagonal energy matrix, and O 
is a orthogonal matrix distributed according to the Haar measure of the orthogonal group. For the 
GOE the distribution of the energies has complicated correlations and a semi circle density, whilc 
they are independently Gaussian distributed for the POE. 

As the level density has no relation to the chaoticity or integrability of the system, we shall unfold 
both spectra to have uniform density with the variance (-E 2 ) normalized to one in the ensemble 
average. As a consequence, the length of the spectrum is 2-^/3 and the average level spacing is 
d = 2^/3/N. One reason to fix the energy scale in this way is, that the variance is easy to control 
even in the weak coupling regime. Note that decoherence depends on the behaviour of the level 
density as a function of energy. Our random matrix model allows also to choose this behaviour 
diffcrently, e.g. as dictated by the classical or quantum Hamiltonian. 

The case of strong coupling both for integrable and chaotic systems was modeled in jïl| with 
two-dimensional anharmonic oscillators, and the systems considered in || might be close to this 
domain. We shall also use equidistant "picket fence" spectra to complete the range of possible 
spectral correlations. The latter is important, because the spectra of low dimensional systems are 
much stiffcr [^6[ [Ï7| than the universal random matrix ensembles would predict. 

The dimensions n, m and N of our Hilbert spaces are chosen finite. If we think of Hamiltonians 
with infinite spectra this truncation is the only reference we make to a choice of the initial wave 
function. 

The entanglement of the two subsystems will be measured in terms of the purity defined as Q 

Z(í) = T ri [Tr 2 (p(í))] 2 = Tr 2 [T ri (p(í))] 2 . (1) 

Here Tri indicates the trace with respect to the first (Latin) index and Tr 2 the one over the second 
(Greek) index. The definition of the purity / is related to the idempoteney defect or linear entropy, 
defined as l-I §. 

We are interested in the time-evolution of an initially pure, i.e. non-entangled state. Therefore 
we construct the initial density matrix /5í Mi í' m '(0) from a "product state" which is pure with respect 
to both pairs of indices, and we have 1(0) = Tri[Tr 2 (p(0))] 2 = 1. Denoting by A the diagonal 
matrix with entries A Q = exp[ií£" Q ] (h is set equal to one), we find in the basis of double indices 
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p(t) = OAO T p(0)OA*O T and by consequence 

I(t) = Tr 1 [Tr 2 (OAO T p(0)OA*O T )Tr 2 (OAO T p(0)OA*O T )] . (2) 

We take the averages involving energies and the averages involving states separately for different 
terms of the sum. In principle we do not need to specify p(0) in morè detail. After performing the 
averages, the result must be independent of the two factor states used to construct p(0). This is 
due to the orthogonal invariance of the total Hamiltonian H used, and the invariance of the purity 
with respect to independent orthogonal transformations in the factor spaces. Hence, without loss 
of generality we may set p(0)u u = 1 with all other matrix elements being zero. With this initial 
condition we obtain the ensemble averaged purity 

= Al,/3;7,<5 B^p-^j (3) 

in terms of the two averages 

= (A Q A 7 A; AJ) = (exp[i t(E a + E 7 - E - E s )]) (4) 

B a ,f3;~t,S = ^ {Oi^^aOii^aOiiJjOj^^Ojv^Oii^On^Oiv^) . (5) 

The averages are connected only because one may force indices to be equal and thus reduce the 
other to a special case; as we shall see below five different terms exist. Note that we omitted the 
average symbols on the purity itself for convenience. 

The obviously relevant time scales are the Heisenberg time d^ 1 = N/(2y/3), and the inverse 
length of the spectrum, which is 1/(2^3). Using these we obtain four different regimes for the time 
evolution, because of the modulus operation inherent in the exponential exp[iiíí]: 

1) Short times, í < 1: Here, perturbation theory can be applied, and we will find the expected 
t 2 dependence with a factor given to leading order by the variance of the energy eigenvalues 
{El) = 1. 

2) First filling of the unit circle at t — n We will find a quadratic minimum for the purity 
with value I mhl = l/n+ 1/m + 0(iV _1 ). 

3) Long times, 2tt/<1 ^> t ^> 7r/v3: In this region, the spectrum has winded many times around 
the unit circle. This acts as a random number generator eliminating correlations. Hence we 
will obtain a result similar to the one for 2) though sub-leading terms may be different. 

4) Poincaré recurrence at t — 2n/d: At this point, a picket fence spectrum will cause exact 
recurrence, while even for a GOE spectrum the recurrence is essentially wiped out. Yet 
for low-dimensional systems with their long-range stiffness ]Ï6| , [r^ | and for models involving 
harmònic oscillators this part may well be important. 

First we shall calculate the energy average A a .p^^- The result does not depend on the vàlues of the 
indices, but only on whether certain indices are equal or not. If two indices of the energies coincide, 
we get either if they have opposite signs or twice the energy if their signs are equal. It may readily 
be seen that five terms are possible (the rightmost equality is vàlid for a random spectrum only): 

5i(í) - (exp[-ií {Ex -E 2 + E 3 - EJ]) = f\t) 
S 2 (t) = (exp[-ií (E 1 -E 2 )])=f 2 (t) 
S 3 (t) = (exp[-ií (2E X -E 2 - Es)]) = f(2t)f 2 (t) (6) 
Si(t) = (exp[-ií 0]) = 1 
S 5 (t) = (exp[-2ií (E, - Et)]) - f 2 (2t) . 
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The first equality in each line determines the special case at hand and the second one gives the 
result for a random spectrum. Thcrc f(t) = sin(v / 3í)/(v / 3í) is the Fourier transform of the level 
dcnsity, which we assumed to be uniform. For GOE spectra the evaluation is morè difficult, but 
some general considerations hold for any kind of spectrum. For long times all terms except S4 go 
to zero. For short times, on the other hand, Si dominates because it has the largest weight. We 
now consider the four time regimes: 

In the short time limit we expand the exponential. Due to the symmetry of the energy distri- 
bution, the linear terms in t vanish while quadratic ones survive. These are of two types. Each 
exponential associated with a given index has a quadratic term, and indices in the linear terms 
of two exponcntials may coincide. This implics that we only need the well-known averages over 
monomials of fourth order in the group elements |ïs|| to obtain 

J(t) ~ 1 - 2{E 2 a ) t 2 [l-(n + m+ 1)/(N + 2)] . (7) 



In the last factor we seem to have al/N correction. Yet if n and to grow as vN the correction is 
of order 1/y/N. If one of the two dimensions is kept constant, the other becomes proportional to 
N, and the second term is of order 1. Terms resulting from correlations of the energies are truly of 
order X/N and were omitted. 

The next time scale is that of the first filling of the unit circle, for which the first minimum of the 
function f 2 (t) is reached. We have a complicated interplay of different terms and it seems that we 
would nccd the average over the orthogonal group B a ^ ]lt s completely. However, this can be avoided 
using the following trick: For uniform density of the spectrum the energy eigenvalues are essentially 
the eigenphases of a circular ensemble. For the case of GOE fluctuations the corresponding enscmblc 
is known as the circular orthogonal ensemble (COE) |ïl"|] , which is the ensemble of unitary symmetric 
N x N matrices S. This ensemble has a unique invariant measure. By identifying the energies with 
eigenphases of a COE the only approximation we make, is that we neglect the correlations that 
exist between the ends of the spectrum for COE. In terms of S we obtain 

/mi„ « (Tii[Ti 2 (S p(0) S*)Tr 2 (S>(0) S*)}) . (8) 

where S^ijv = J2 a 0,· í4) aexp[iS a (27r/r)]0 J · 1/)a . The ensemble average originally given as one over 
states and spectra is thus given in terms of averages over four symmetric unitary COE matrices, 
two of which are complex conjugate. Such averages are calculated in [Ï9[ | and we obtain 

(n + m)N 2 + [3(n + to) + 2] N - 2(n + to - 1) 
min ~ N(N + l)(N + 3) ■ { ) 

As we shall see below this is slightly lower than the long time limit, while for the integrable (POE) 
case the long time limit and the value J m i n coincide. We shall calculate the long time limit next. 

For t ^> 1, the process of stretching and taking modulo 2ir is a reasonably efficient randomizer for 
a fluctuating set of numbers with correlations such as a GOE spectrum. Therefore the eigenphases 
on this time scale are random both for the GOE and the POE. Thus only the fourth term survives, 
where the indices of energies in conjugate terms coincide. The energy dependence, and therefore 
the time dependence, drops out and we are left with averages over the orthogonal group. Only 
two- vector terms i.e. averages over elements from two rows of the matrix oceur. These have been 
calculated poj , and we find 

_ (n + m)N 3 + 3[4(n + m) + 3)N 2 + [35 (n + to) + 57}N + 48 
°° ~ (N + l)(N + 2)(N + 4)(N + 6) ' ( } 

For the POE this result holds equally at time t = 7r/v3, which we have discussed above for the 
GOE, though it will oscillate for larger times. The value for the GOE lies slightly below the one 
for POE. What we see is a weak signature of the correlation hole, characteristic of chaotic systems. 

When evaluating the large N limit of these expressions, we have to take into account the n and 
to always oceur in the form (n + m). As N = nm goes to infinity, n + m can behave as any power 
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Figure 1: The average purity as a function of time for m = n = 4. The ensembles considered are: 
POE (long dashed line), GOE (sòlid line), and picket fence (short dashed linc). The value Ioo is 
denoted by a thin dotted horizontal line; the value I m i n at the first minimum by a thin upright 
cross. We clip a small interval of the ordinate, to see the difference between GOE and POE in more 
detail. To observe the behaviour of the purity at short times, see Fig. 2. 



N q , with 1/2 < q < 1. The two extremes are realized if n and m increase simultaneously, such that 
the ratio n/m remains constant, or if one is fixed such that the other becomes proportional to N. 
Keeping this in mind we can still expand both expressions and find that the first minimum for the 
GOE case will always be slightly lower, but the difference will diminish as Í/N 2 , while purity itself 
will either diminish as 1/y/N or reach a constant value í/n or í/m. We can therefore conclude 
that the effects of spectral correlations are quite insignificant except for the smallest dimensions in 
both subsystems. 

For times of the order of the Heisenberg time we expect a very different behaviour. More 
precisely, \ï t = 2n/d = Nn/y/3, we have exact revival in the case of the a picket fence spectrum, 
while there will be no particular signature for a random spectrum. In the case of GOE fluctuations 
one should consider the width of the kth neighbour spacing distribution. It is known to increase 



logarithmically. For k = 1 it has a width of w 1.25<i and for k = 8 its width already is « 1.85<i |21 



Due to the increasing broadness, we will see no Poincaré revival in the case of a GOE spectrum 
either. We should though note two facts: First we will find an additional partial revival at half the 
time mentioned for a picket fence spectrum, because of the terms S^í) and Ss(t), which contain 
eigen-energies with a factor two. Second, and more important, the long-range stiffness of spectra 
in low-dimensional systems |Ï6j implies a saturation of the width of the fcth neighbour spacing 
distribution and could therefore lead to recurrence effects. 



4 Numerical results 

First we consider a strong coupling situation with small dimensions, where all the effects we predict 
are most notable. For this purpose we choose the case n = m = 4 and show the corresponding 
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Figure 2: The average purity as in figure 1, but for the GOE case only (thick sòlid line). The 
average purity plus and minus the Standard deviation (dotted lines), and the purity for the first 
mcmber from the ensemble (dashed line). 



time evolution of purity in Fig. 1. The results are displayed for spectra with random and GOE 
like fluctuations as well as picket fence spectra. The ensemble size is 2 * 10 6 . As expected, at the 
beginning all curves are equal, and therefore we choose a scale where the quadratic dependence 
at the origin is not visible. The theoretical predictions are well fulfillcd. The first minimum 
occurs at t = tt/VÏ for all spectral ensembles. Its depth for GOE correlations coincides with the 
COE result indicated by a cross, while both the GOE and the POE case take the asymptotic 
value given by a dotted line. In the POE case, the oscillations are essentially due to the form of 
f{t) = sin(v / 3í)/(v / 3í) (cf. eq. (||)). However, the GOE and picket fence cases follow this behaviour 
in a similar fashion. It is essentially a "diffraction" effect due to the sharp cutoff of the level density. 
E. g. for a semi-circle density the effect is attenuated, and for a Gaussian density it disappears. 
Therefore the experimental significance of these oscillations is limited to special situations, but the 
ràpid decay on the time scale indicated does not depend on this fact. The first minimum is lower 
for GOE like fluctuations and for picket fence spectra than for random spectra, though the effect 
is only a few percent. The results coincide with our theoretical predictions both at the minimum 
and in the asymptotic region. This is also true for the recurrences, which are only seen for picket 
fence spectra both at t = and more weakly at t = Nn/ (2v3)- 

If we neglect oscillations, the rise of purity after the first minimum follows roughly that of 
the Fourier transform of the two-point function with appropriate scaling. This is not surprising, 
because we may expect that a cluster expansion of the correlations relevant for the difference from 
the random case is dominated by the two-point function. Yet it is important to note that the 
diffcrences between the POE and GOE case lie within a few percents for such small dimensions. 
If n or m is increased any further, the differences become practically invisible. In || we show a 
figure for n = m = 10, which confirms this fact. We conclude that decoherence of random states is 
inscnsitive to chaos or integrability, except for very small systems. 

Concerning noticeable differences between the integrable and chaotic cases, the news get even 
worse, if we study the fluctuations of the purity as a function of time. Fig. 2 shows the GOE 
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case for n = m = 4 as above but including the Standard deviation of the purity plotted as a band 
around its average. For illustration we have also included the first mcmber of the ensemble into the 
figurc. Thc Standard deviation is big, and it is fair to say, that any medium and long time features 
shown by a particular initial state in a chaotic system are most likely fluctuations. The picture for 
a random spectrum is essentially the same and thus the same conclusion holds for a random wave 
function in an integrable system. 

Calculations with larger dimensions show that thc fluctuations of the purity diminish, but not 
as fast as the difference resulting from spectral correlations. This is shown to diminish with N~ 2 . 

Finally let us discuss some numerical results for weak coupling: We construct our ensemblcs in 
the eigenbasis of . Then fcW and are diagonal with elements taken from the spectrum 
of the GOE or the POE. In both cases, thc spectra arc unfoldcd to uniform density. The matrix 
elements of l/( 1,2 ) are Gaussian distributed random variables. As the diagonal elements are set 
to zero, we have: ((V^i' 2 ^) 2 ) = 1 — 5 a p. It is assumcd, that thc diagonal elements are effectively 
absorbed into the Hamiltonians and h^ 2 \ 

The resulting level density of the total system becomes quite complicated. However, it is still a 
simple task to obtain its variance: 

(E 2 a ) = ±TrH 2 = (e 2 ) + (e 2 ) + \ 2 {N 1) (11) 

where (ef ) and (e 2 ) denote the variances of the spectra of the subsystems, which are equal to onc 
half. In what follows, the energy scale is again chosen such that (i? 2 ) = 1. 

While the randomness of the initial states was irrelevant in the previous examples of strong 
coupling (due to the orthogonal invarianec of thc total Hamiltonian) , here it makes an important 
difference. In distinction to eigenstates of or localized wave packets, we may expect that 

random initial states decohere faster. The reason is most easily understood looking at the LDOS. 
Thcrc a random initial state overlaps typically with many more eigenstates of the total system, as 
localized wave packets or eigenstates of H^°\ In the course of time this leads to faster decoherence. 

Again small dimensions provide us with some surprises: Fig. 3 shows results for the dimensions 
n = m — 4. The average purity is displayed for three different combinations of spectra used in ftW 
and h {2 \ namcly for GOE-GOE, GOE-POE, and POE-POE. The coupling parameter is A = 0.03. 
We find that for short and intermediate times the three curves have the same behaviour, but the 
asymptotic vàlues diffcr: In contrast to the intuitive picture, that decoherence should be strongest 
in the chaotic case, the asymptotic value of the purity is lowest for the POE-POE case, which 
corresponds to two integrable subsystems. We get an intermediate value for the GOE-POE case, 
and the highest value is obtained for the GOE-GOE case, corresponding to two chaotic subsystems. 

Yet again the significance of this difference is limited. This becomes clear by means of Fig. 4. 
Thcrc we show again the Standard deviation of thc purity in a similar fashion as in Fig. 2. This 
is done for the GOE-GOE case from Fig. 3. In addition the purity curve for a single Hamiltonian 
from the ensemble is plotted. As can be seen in the figure, the size of the fluctuations of this curve 
is very well described by the Standard deviation. 

Wc sec that the band is extremely wide and much bigger than thc difference in the purity 
between the different combinations of spectral statistics. Note that the corresponding figure for the 
other two cases from Fig. 3 (not shown) yicld thc same conclusion. 

In Fig. 5 we show again the average purity for the three different combinations of spectral 
statistics in the subsystems, but here the dimensions are n = m = 10, and the coupling parameter 
is A = 0.01. The smaller coupling parameter (as compared to the previous case, where n = m = 4) is 
meant to compensate approximately the increase of the total dimension TV. Note that the variance 
of the spectrum of \V^ 1,2 \ which may serve as a measure for the strength of the perturbation, 
increases linear ly with N. 

Thc results shown in Fig. 5, are very similar to the previous ones. However, the differences 
between the three curves at large times has diminished. Note that an increase of A would further 
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Figure 3: The average purity as a function of timc for n = m = 4, in the weak coupling case 
(A = 0.03). It is plotted for three different combinations of ensembles for M 1 ) and h^: GOE-GOE 
(sòlid linc), GOE-POE (dashed line), and POE-POE (dotted linc). 




t 



Figure 4: The average purity as in Fig. 3, but for the GOE-GOE case only (thick sòlid line). The 
average purity plus and minus the Standard deviation (dotted lines), and the purity for the first 
member from the ensemble (dashed line). 
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Figure 5: The average purity as a function of time for n = m = 10 in the weak coupling case 
(A = 0.01). It is plotted for three different combinations of ensembles for hP^ and h^: GOE-GOE 
(sòlid linc), GOE-POE (dashed line), and POE-POE (dotted line). 



reduce these differences. The Standard deviations of the purity curves also diminish (not shown), 
but still exceed the differences in the asymptotic average vàlues. Hence, there is probably no chance 
to see an effect of the different spectral statistics of the subsystems without doing averages over 
rather large samples. Nevertheless a detailed analysis of the different decays in weak coupling is an 
interesting open question. 



5 Conclusions 

The influcncc of integrability versus chaos of the dynamics of a system on the entanglement of a 
random product wave function was considered. The type of dynamics was expressed in terms of 
spectral statistics, as this is the only invariant property of quantum dynamics, which is sensitivc 
to the above distinction. We conclude that the effect of the type of dynamics on random states 
is small, and vanishes rapidly with the dimension of the Hilbert spaces involved. For quantum 
information processing, where random states are typical, this means that the distinction between 
chaos and integrability is not relevant to the rate of entanglement of two subsystems and thus for 
the decoherence in the central system. 

Some effects have been seen for small systems, and these are interesting in themselves. In 
particular we find very large fluctuations for individual purity decays. One may argue that the 
fluctuations seen for initial packets in the chaotic area of phase space in [|| , could be interpreted as 
fluctuations, but clearly a large statistical calculation would be necessary. 

A comparison with the work of Prosen and Znidaric |^|, [| is more difficult. As the Hamiltonian 
is disturbed we are looking at non-linear effects even in quantum mechanics. After all this is what 
the fidelity concept was designed for. In this context our result is important because, the reduction 
in the loss of fidelity, that is obtained in ^ || for chaotic dynamics, does not imply increased 
decoherence in the forward time evolution scheme with a fixed Hamiltonian, as one might have 
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inferred from the common concept, that chaos enhances decoherence. 
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